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Jj 要 : 针对 间 歌 过 程 的 非 线性 和 动态 性 ， ee (GLRGMM) 算 法 。 首 先 引 入 令 域 


保持 诺 入 算法 提取 局 部 流 形 结构 ， 通 过 寻求 一 种 低 维 投影 对 非 线 性 进行 全 局 结构 保持 ， 同 时 最 大 限度 地 保留 局 部 
流 形 特征 ; 然后 通过 对 高 斯 混合 模型 引入 正则 项 来 在 线 est 获取 非 线 性 数据 流 形 结 构 ， 解 决 数据 动态 
性 问题 ; 最 后 集成 全 局 一 局 部 监控 指标 实现 在 线 监 控 。 通 过 青 者 素 发 酵 过 程 进 行 了 验证 ,结果 表明 所 提 算 法 比 DPCA、 


GLNPE 具有 更 好 的 在 线 监 控 效 果 。 
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Online monitoring strategy for batch process based on GLRGMM 
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: Abstract: Aiming at nonlinear and dynamic characteristics of batch process, this paper proposed a global-local regularization 

© Gaussian mixture model (GLRGMM) algorithm. At first, this algorithm introduced Neighborhood Preserving Embedding 

Q algorithm to extract the local manifold structure. It preserved global structure of nonlinear process by seeking a 
low-dimensional projection while preserved local manifold features to the maximum extent at the same time. Then it 
introduced the regularization term into GMM to online monitor and update Gaussian model. It obtained the manifold structure 
of nonlinear data and solved the data dynamic characteristic problem simultaneously. Finally, it used the integration of 
global-local monitoring indicators to effectively achieve online monitoring. The results verified by the penicillin fermentation 
process show that the proposed algorithm has better online monitoring effect than DPCA and GLNPE algorithms. 
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间歇 过 程 中 的 非 线 性 和 动态 性 问题 , KPCAU!, DPCA 和 DPLSP! 

应 运 而 生 并 得 到 了 成 功 应 用 。Lee 等 人 名 提 出 了 多 向 核 主 成 分 

间 鞭 过 程 由 于 其 操作 性 强 和 灵活 度 高 等 优点 成 为 一 种 重要 分 析 算 法 (multi-way kernel PCA, MKPCA); Peng 55 A! fe T-PLS 

的 工业 生产 方式 ， 特 别 是 在 精细 化 工 、 生 物 制药 和 食品 饮料 等 ”的 基础 之 上 提出 了 全 影 核 偏 最 小 二 乘 (total kernel partial leasts 
小 批量 、 高 附加 值 产品 的 生产 和 制备 当中 得 到 广泛 应 用 上。 强 quares，T-KPLS) 算 法 。 这 些 算法 在 一 定 程度 上 能 有 效 解 决 过 程 
非 线 性 、 动 态 性 是 间 葡 过 程 的 固有 特征 ， 是 实际 生产 过 程 和 在 变量 之 间 的 非 线 性 特性 ， 但 仍 存在 计算 量 大 和 核 参 数 选 择 盲目 
线 监控 中 面临 的 问题 ， 而 寻求 一 种 有 效 的 解决 途径 对 降低 经 济 性 问题 ， 并 且 忽 略 了 时 序 相关 性 ， 没 有 考虑 动态 性 问题 ， 难 以 
损失 和 避免 人 员 伤 亡 具有 重要 意义 中 。 保证 过 程 监控 的 准确 性 。Fan 等 人 外 提出 了 基于 核 函 数 的 动态 
随 着 工业 过 程控 制 系统 的 发 展 ， 大 量 过 程 数据 被 获取 并 得 独立 元 分 析 KDICA 算法 ; Lu 等 人 外 提出 了 双向 动态 主 成 分 分 
到 有 效 利 用 ， 其 中 以 PCA、ICA、PLS 为 主要 代表 的 基于 多 元 析 方 法 。 这 些 方法 都 是 通过 历史 数据 扩展 当前 数据 来 消除 时 序 
统计 过 程 监控 (multivariate statistical process monitoring, MSPM) 相关 性 ， 并 不 能 很 好 地 消除 过 程 数 据 的 动态 性 。 为 了 同时 解决 
方法 在 工业 生产 过 程 中 得 到 了 迅速 发 展 和 广泛 应 用 。 为 了 解决 过 程 数据 中 的 非 线 性 和 动态 性 ，Wang 等 人 00 提 出 了 双向 动态 
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核 主 成 分 分 析 法 (dynamic kernel PCA, DKPCA); Zhang 等 人 0 
提出 了 递归 核 PCA 方法 ; Jia 等 人 0 提出 了 基于 批 动态 核 PCA min > ， 
方法 。 这 些 方法 虽然 能 够 同时 监测 过 程 数 据 非 线性 和 动态 特性 ， 


X 2x 
j 


; Dw; =1,7=1,2,.…,N (1) 
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但 没有 考虑 过 程 数据 的 全 局 和 局 部 结构 信息 ， 仍 存在 主 元 个 数 映射 变换 矩阵 A 通过 式 OO 最 优 解 得 到 。 
多 、 故 障 误 报 率 较 高 等 问题 。 : 

近年 来 ， 流 形 算法 由 于 其 优越 的 局 部 流 形 结构 特征 提取 和 E Y Yn [ conie o 
高 维 数据 处 理 能 力 而 越 来 越 多 地 应 用 于 过 程 监控 领域 .Hu 等 人 
应 用 局 部 保持 投影 算法 (locality preserving projection, LPP)U?!fü NPE 虽然 能 很 好 地 保持 局 部 流 形 结构 ， 并 能 揭示 隐藏 在 高 


AU B I FFA SE iE (neighborhood preserving embedding, NPE)'^! 维 数据 中 的 非 线 性 ， 但 忽略 了 全 局 结构 ， 并 不 能 更 好 地 保持 数 
对 间歇 过 程 进行 了 有 效 监 控 ;， Zhang 等 人 05 提 出 了 一 种 基于 全 据 结构 信息 。 因 此 ， 本 文 首先 提出 了 全 局 -局 部 邻 域 保持 嵌入 
局 一 局 部 结构 分 析 (global-local structure analysis, GLSA) 算 法 ， (global-local neighborhood preserving embedding,GLNPE) 算法 。 
保持 了 全 局 结构 信息 的 同时 获取 了 局 部 流 形 结构 :Ya09 根 据 1.2. 全 局 一 局 部 邻 域 保持 戏 入 (GLNPE) 算 法 

GLSA 提出 了 局 部 一 全 局 PCA (local and global PCA, LGPCA) 1.2.1 全 局 结构 保持 

算法 ; Tong 等 人 0 提出 了 多 流 形 投影 (multi-manifold 
projection, MMP) 算 法 , 结合 了 全 局 和 局 部 方差 信息 用 于 过 程 监 


对 于 R ZERE X [exa] ， 寻 找 投影 矩阵 


veo 控 。 这 些 算法 虽然 提取 了 数据 的 全 局 和 局 部 信息 ， 但 并 没有 有 A ， 通 过 投影 映射 了 = AT X 将 其 投影 到 低 维 特征 空间 ， 根 据 
之 效 解决 数据 的 非 线 性 和 动态 性 问题 。 由 于 高 斯 混合 模型 方差 最 大 方向 来 最 大 程度 地 保留 原始 数据 全 局 结构 信息 ， 其 
(Gaussian mixture model, GMM) 可 以 更 好 地 处 理 复 杂 过 程 数 据 标 函 数 如 式 (3) 所 示 。 
而 被 用 于 估计 复杂 数据 分 布 ， 并 成 功 应 用 于 非 线性 、 动 态 特性 Ls (A) siy beni -XJ A- A'CA a 
B3 TALEE lS- nm 
为 了 解决 间 葡 过 程 的 非 线 性 和 动态 特性 ， 并 且 有 效 提取 过 n 


1 
其 中 : z-lyx, C=-2(% x)x-x). 
n 


程 数据 的 全 局 一 局 部 信息 ， 本 文 提 出 了 全 局 一 局 部 正则 化 高 斯 m 

混合 模型 (global-local regularization Gaussian mixture model, 1.2.2 局 部 结构 保持 

GLRGMM) 算 法 。 由 于 NPE 可 以 更 好 地 获取 高 维 数据 的 局 部 流 局 部 结构 保持 是 通过 每 个 数据 点 及 其 近邻 点 进行 重 构 来 控 
T 


形 特征 ， 通 过 对 NPE 进行 改进 求 取 一 个 新 的 低 维 投影 ,获取 高 。 掘 出 数据 局 部 结构 信息 ， 其 保留 局 部 结构 的 目标 函数 如 式 (4) 


维 空间 全 局 -局 部 流 形 结构 , 实现 非 线 性 过 程 数 据 降 维和 特征 提 所 示 。 
取 ; 然后 对 GMM 进行 动态 流 形 正 则 化 ， 自 动 更 新 获取 当前 数 ; 2 
据 ; 最 后 通过 NLLP (negative log likelihood probability) fll MD J (A)= min 》， M -Ywy =min7 (I -W) (1-W)Y 7 
(mahalanobis distance) 统 计 全 局 、 局 部 概率 生成 综合 监控 指标 u j " " s P x 
GLP (global local probability)， 实 现 对 间歇 过 程 的 在 线 监控 。 三 Xx(r-W) (I-W)X 8-352 XMX A 
1 ”全 局 一 局 部 邻 域 保 持 嵌 入 算法 其 中 ，W 为 权重 矩阵 ，M =(1-W) (1-W). 
1.1 邻 域 保持 诅 入 算法 (NPE) 1.2.3 全 局 一 局 部 邻 域 保持 说 入 算法 (GLNPE) 
NPE 算法 能 够 挖掘 出 高 维 空间 中 隐藏 的 低 维 流 形 结构 ， 能 为 了 更 加 充分 地 提取 过 程 数 据 特征 ， 需 要 找到 一 种 投影 使 


保持 原始 空间 的 局 部 流 形 结构 。 其 目的 是 将 R?” 空间 的 高 维 数 。” 得 在 低 维 空间 同时 保留 全 局 和 局 部 信息 ， 即 在 式 (3) 最 大 的 同 
时 能 保证 式 (4) 最 小 。 为 了 解决 这 个 复杂 优化 问题 ， 本 文 提出 
了 新 的 目标 函数 : 


据 映 射 到 一 个 相对 低 维 的 特征 空间 了 = ye Yon Xa] o 


{ylyeR(d<D) ， 同 时 求 出 相应 的 变换 矩阵 | Ju (A) 
J (A) = min 
NJ siobal (A)+ (1-7) Joca (A) (5) 
A - (a2, 7,24] . : AT'XMX'A 


OUT gA'CA* (1-5) A XMX' A 


首先 根据 大 近邻 法 确定 样本 点 大 的 大 个 近邻 点 并 组 成 局 部 


其 中 : 0<7 <1 为 可 调 平衡 参数 ， 本 文 取 7 =0.85 。 式 (5) 
最 小 化 问题 可 以 通过 拉 普 拉 斯 方法 求 得 投影 矩阵 A : 


RM (x), ， 再 通过 近邻 点 来 计算 权重 系数 矩阵 凤 ， 权 重 系 


数 矩 阵 可 通过 求解 式 〈1) 的 最 优 解 得 到 。 


(7C+(1-7)M)4A=44IXMXTA4 (6) 
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此 可 得 低 维特 征 空间 数据 : n 
LAT Y s,»(c.|x;) 
Y-AX (7) s) (s-1) j=1 
"H3 " Pr 5B Pa B: Hd ST P NÀ = P (C. |x;) E p (c, |X;)- n (13) 
通过 投影 映射 A 将 高 维 数据 投影 到 低 维特 征 空间 的 同时 Ys, 
最 大 程度 地 保留 了 全 局 一 局 部 结构 信息 ， 更 有 助 于 对 非 线性 数 HC 
据 特 征 的 提取 。 M-step 
2 ”全 局 一 局 部 正则 化 高 斯 混合 模型 (GLRGMM) 算 法 » p" (c, |x;) T 
X^? i-i 
n 
[Ett dae fts e MR D-(X,]i2 L2, n), Xn 
n (s) 
为 变量 个 数 ， 每 个 变量 的 维 数 为 加 ， 则 样本 的 高 斯 混合 模型 本 | 
Fen py (15) 
. 之 (Cj 
P(X)= Do.8 (XB) (8) 
k=1 n T 
(s) E (s+1) (s+1) 
Rep. Ko, AZ 分 别 表示 高 斯 成 分 数 ， 第 个 高 斯 成 分 的 yen 2p (c. x. (x, D "(x 一 从 ) a 
k n 


权重 、 均 值 和 方差， 其 目标 函数 可 以 写成 


J (8) = max] log p(X) |= max jal To(s e) (9) B 


erp. oC, ui x90 yg dese LI rp B Ic 


p* (c, |X;) 


通过 GLNPE 提取 特征 


空间 中 的 两 个 数据 点 羡 ; 和 X; ， 大 


余弦 相似 性 甜 阵 8 = cos Z(X,, X;) . 成 分 的 权重 、 均 值 和 方差 。 
GLRGMM 通过 GLNPE 对 非 线 性 数据 的 全 局 一 局 部 结构 
定义 正则 项 R, : RM m 
特征 进行 了 更 加 充分 地 提取 ， 更 符合 实际 间歇 过 程 数 据 特性 。 
1 其 详细 的 算法 流程 如 图 1 所 示 。 
im DIG (c6. |X.)- plc IX, ) s, (REA 4E Y 
2 P 1 (样本 集 x ) 
-Y»c. Ix.) 1, DUC Ix)p(G.[X;)S; ao 数据 标准 化 
E fe Life -fo Sf. GLNPERFT-RER 
- fc He Y Y 
全 局 结构 保持 | | 局 部 结构 保持 | 
其 中 ，; 为 单位 矩阵 ， L-I-S; p(C.|X;) 为 X; 属 于 第 k 
正 数据 Y 
AN Es HE pi Ly ft 日 
个 高 斯 成 分 的 后 验 概率 ， RECITE 
ir Y 
fi =| p(C.|X),…,p (ClX, )| (11) K-means 初 始 化 参数 0 
通过 最 小 化 Ri 确定 后 验 概率 , 使 得 不 同 高 斯 成 分 之 间 平 滑 | 在线 更 新 后 验 概率 四 (GX) 
过 渡 ， 去 掉 后 验 概率 最 小 的 高 斯 成 分 ， 更 新 当前 高 斯 成 分 ， 实 "ES ap? WE 
现 动 态 过 程 平滑 过 渡 和 数据 实时 更 新 ,引入 一 个 包含 正则 项 R 
EE SEPTA NDS 
的 新 目标 函数 : J 
p ui E N 
1,.(6)- yp (C. |x, Jlog (os e (X, lu. X .)-48 (12) - AE 0) I 973 3 
i-l k=l = 
Jr. 4 为 正则 化 参数 。 采 用 K-means 算法 得 到 初始 参数 -— — 
、 i ue 返回 GLRGMM 人 参数 0 
K,,0,, 14,2. 通过 EM 和 帮 代 自动 更 新 确定 高 斯 成 分 参数 
图 1 GLRGMM 算法 流程 
(K,o,u,X). Fig.l Flow chart of GLRGMM algorithm 
E-step: 3 #F GLRGMM 的 监控 策略 
通过 最 小 化 R, 确定 后 验 概率 : 
通过 最 小 化 R, 确定 后 验 概率 3.1 间歇 过 nesem 
[R] Sk xc fe Sk mp Nm NU f su mu e dk 7g — EU: 
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X(xJxK), ， 其 中 ，7 为 批 次 个 数 ，v 为 变量 个 数 ， 天 为 


采样 时 间 。 在 实际 的 生产 过 程 中 ， 由 于 不 同 批 次 的 采样 时 间 不 
同等 原因 而 不 可 能 做 到 完全 重复 同 批 次 的 生产 ， 这 就 不 可 避免 
地 出 现 数据 不 等 长 问题 。 本 文采 用 最 短 长 度 法 20， 以 最 短 的 一 
个 批 次 数据 作为 标准 而 将 其 余 批 次 数据 进行 硬性 截取 ， 使 得 所 
有 批 次 数据 都 具有 相等 长 度 。 

对 间歇 过 程 的 三 维 数据 进行 建 模 ， 需 要 把 已 处 理 成 相等 长 
度 的 数据 展开 成 与 连续 过 程 相 类 似 的 二 维 数据 。 本 文 结合 了 批 
次 展开 可 以 去 除非 线性 、 变 量 展开 可 以 消除 数据 预 估 问 题 的 优 


点 ， 展 开 方 式 如 图 2 所 示 。 
K 
Timeo AAA Qd 
ox di k=1 
| 1x2 
1 fe—— ZEE 21 —— 
Batch) I : 
k=d 
T : 
Variable(/) KI 
i J £2J kJ 
1 
k=1| k=2| k-d 
I 


图 2 间 敬 过 程 数据 展开 方式 
Fig.2 Scheme of batch process data unfolding 


首先 , 将 三 维 数据 X (1xJx) 按 批 次 方向 展开 成 二 维 数 
ik X (1x JK) 并 进行 标准 化 处 理 : 

- 1 ~ 
Xj g I 2 X jk (17) 

i-l 
ou = (18) 

E XPik T X 

X ik 一 (19) 


得 到 的 数据 在 一 定 程度 上 可 以 消除 非 线 性 和 动态 特性 ， 然 


i 


置 为 三 维 数 据 


后 将 标准 化 的 二 维 


数据 X xJK) 8 


展开 成 二 维 数据 X (KI xJ) 并 进 


X (IxJx K), 再 按 变 量 方 


a 


行 标准 化 ,这 种 展开 方式 可 以 消除 在 线 监 控 时 的 数据 预 估 问 题 ， 
提高 了 实时 监控 的 精确 率 ， 降 低 了 误 报 率 。 
3.2 GLP(global local probability) 监 控 指 标 
为 了 统计 全 局 概率 指标 ， 首 先 引 入 负 对 数 似 然 概 率 
Cnegative log likelihood probability, NLLP). 


NLLP(X.,) - -log p(C, |X;) (20) 


全 局 概率 指标 通过 


比较 测试 数据 与 训练 数据 的 NLLP 得 到 


P, (X,,) - Pr[NLLP(X 


8 


)SNLLP(X,,,)) (2D 
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引入 服从 好 分 布 的 马 氏 距离 (Mahalanobis distance, MD): 


MD(X,)-(X,-u)(Z,) (X; -4y (22) 


部 概率 指标 : 


pO (x, )= Pr[MD((X 


l 


定义 局 


ClX € C,) 


train ? train 


(23) 
< MD((X,,.C, )|X,., e C, 


test? 


局 部 概率 生成 GLP(global local probability) 


监控 指标 : 
OLP =Ý P, (CXu) E" (X...) 24) 
33 ”离线 建 模 
a) 在 正常 工 况 下 采集 若干 批 次 数据 匀 (TxJ x 天) 作为 训 


练 样本 ， 采 用 3.1 节 的 方法 进行 数据 预 处 理 ， 最 终 得 到 已 标准 


化 处 理 的 二 维 数据 XX(KIxJ); 


局 部 信息 保持 ， 由 式 (6) (7) 


b) 通过 GLNPE 进行 全 局 


得 到 低 维特 征 空间 投影 数据 了 ; 
c) 利用 GLRGMM 对 数据 了 进行 正则 化 ， 用 EM 算法 确 
定 参数 0 ; 


d) 通 过 式 (20)(22 ) 分 别 计算 NLLP(X pain) > MD(X pain ) 。 


3.4 在 线 监 控 
a) 将 测试 数据 按 3.1 节 的 方法 进行 数据 预 处 理 。 
b) 通过 GLNPE 算法 而 得 到 新 的 测试 数据 的 投影 数据 了 。 
c) 通过 GLRGMM 对 新 的 投影 数据 了 进行 正则 化 , 用 EM 
算法 确定 参数 0 。 


d) 通过 式 (20)(22) 分 别 计算 NLLP (X pu) , MD(X,,) 。 
e) 通过 式 (21) (23) 分 别 计算 全 局 、 局 部 概率 指标 。 


过 式 (24) 计算 GLP 监控 指标 ， 若 GLP>0.99， 则 出 
否则 ， 重 复 步 骤 b) ~f) 直到 出 现 故 障 报警 或 生产 过 


现 故 障 ; 
程 结束 。 


4 案例 研究 


青 才 素 作为 一 种 重要 抗生素 而 被 临床 广泛 应 用 ， 其 生产 制 
备 过 程 是 典型 的 间歇 过 程 ， 有 具有 强 非 线性 、 动 态 性 等 特性 。 本 
文 基 于 美国 Ilinois 州立 理工 学 院 开 发 的 Pensim 2.0 标准 青霉素 
仿真 平台 呈 9 产生 出 间歇 过 程 数据 ， 它 可 产生 出 不 同 初始 条 件 和 
不 同 工 况 下 青霉素 发 酵 过 程 中 各 变量 每 个 时 刻 的 数据 用 于 分 析 
研究 。 

本 文 将 每 批 次 的 反应 时 间 设 定 为 400 h， 采 样 时 间 设 为 0.5 
h, 在 正常 范围 设置 不 同 初始 条 件 , 不 引入 故障 的 情况 下 共产 生 
35 个 批 次 正常 工 况 下 数据 ， 从 产生 的 18 个 变量 数据 中 选择 其 
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中 10 个 过 程 变量 作为 监控 变量 〈 表 1)， 构 成 三 维 数据 


X(35x10x800) 作 为 训练 样本 。 


dl 过程 变量 


Table 1 Process variables 
变量 号 变量 名 称 
1 通风 速率 (L/h) 


2 搅拌 功率 (w) 

3 补 料 温度 CK) 

4 溶解 氧 浓 度 (mmol/L) 

5 反应 器 体积 CL) 

6 排 气 二 氧化 碳 浓 度 (mmoyL ) 
7 


pH 值 
8 AREE K) 
9 TER K) 
10 冷水 流 加 速率 (L/h) 


Pensim 2.0 提供 了 三 种 故障 类 型 (通风 速率 、 搅 拌 功 率 、 
底 物 流 加 速率 )。 本 文 引入 故障 类 型 1 ( 即 变量 1 通风 速率 ) f 
故障 类 型 2( 即 变量 2 搅拌 功率 )， 在 时 间 200~300 h( 即 采样 点 
400~600) 加 入 15% 的 阶 跃 信号 作为 故障 信号 ， 所 产生 的 数据 件 
为 故障 测试 数据 以 供 在 线 监控 。 通 过 GLRGMM 对 训练 数据 建 
模 得 到 高 斯 成 分 如 图 3 所 示 。 分 别 用 DPCA. GLNPE f 


T 


nt 


H 


GLRGMM 算法 进行 监控 ， 对 于 DPCA 和 GLNPE 给 出 故障 1 


的 监控 图 ， 如 图 4-7 所 示 ; 其 余 详细 监控 结果 见 表 2、3。 


GLRGMM 对 故障 1 和 2 的 监控 结果 如 图 8、9 所 示 。 


6 0 O 
训练 数据 
正常 测试 数据 
故障 数据 
三 2 
z 
[o] 
£5 
o0 
-2 
E a 0 1 2 3 
GLRGMM 
图 3 “训练 数据 建 模 结果 
Fig.3 Modeling results of training data 
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A] 


4 DPCA 对 故障 1 的 7? 监控 图 


Fig.4 T? monitoring diagram of DPCA for fault 1 
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图 5 DPCA 对 故障 1 的 SPE 监控 图 


Fig.5 SPE monitoring diagram of DPCA for fault 1 


120: 
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图 6 GLNPE 对 故障 1 的 7? 监控 图 


Fig.6 T? monitoring diagram of GLNPE for fault 1 


GLNPE-SPE 
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图 7 GLNPE 对 故障 1 的 SPE 监控 图 


Fig.7 SPE monitoring diagram of GLNPE for fault 1 
表 2 不 同方 法 对 故障 1 的 监控 统计 /% 


Table 2 Monitoring statistics of different methods for fault 1 /% 


方法 DPCA DPCA GLNPE  GLNPE GLRGMM 
(统计 量 ) CT?) CSPE) CT?) (SPE) (GLP) 
误 报 率 26.36 10.7 3.02 9.18 0.5 
漏 报 率 19.4 16.92 3.98 11.44 0 
检测 率 80.6 83.08 96.02 88.56 100 
1 
0.8 | 
0.6 | | 
可 | TN | | 
i | | | | | 
T 中 | | | 
0.2 | | Nl | 
J| | | | | | | | 
MRNA E, Li 
% 200 400 600 800 
Samples 


图 8 GLRGMM 对 故障 1 的 监控 图 


Fig.8 Monitoring diagram of GLRGMM for fault 1 


录用 定稿 赵 小 强 ， 等 : 一 种 基于 GLRGMM 的 间 歌 过 程 在 线 监 控 策略 


参考 文献 : 


"|l [1] RSE, X369], FA. 一 种 不 等 长 的 多 模 态 间 歌 过 程 故障 检测 方法 
a 981 | UJ]. 化 工学 报 , 2016, 67 (7): 2916-2924. (Guo Jinyu, Yuan Tangming, Li 
ii 94r | | | Yuan. Fault detection method for uneven-length multimode batch processes 
0.2 | | [J]. CIESC Journal, 2016, 67 (7): 2916-2924. ) 
o : | [2] £4, Sd, mA. 等 . 基于 角 结 构 统 计量 的 MKECA 间 歌 过 程 故 
0 200 400 600 800 
Samples 障 监 测 DJ]. 仪器 仪表 学 报 , 2017, 38 (1): 174-180. (Wang Pu, Li Chunlei, 
图 9 GLRGMM 对 故障 2 的 监控 图 Gao Xuejin, et al. Research on batch process monitoring method based on 
Fig.9 Monitoring diagram of GLRGMM for fault 2 multi-way kernel entropy component analysis and angular structure 
363 不 同方 法 对 故障 2 的 监控 统计 /% statistic [J]. Chinese Journal of Scientific Instrument, 2017, 38 (1): 
Table 3 Monitoring statistics of different methods for fault 2 — /96 174-180. ) 
方法 DPCA DPCA GLNPE  GLNPE GLRGMM [3] 起 小 强 , 王涛 . 基于 TGNPE 算法 的 间 歌 过 程 故障 诊断 [J] 化 工学 报 ， 
〈 统 计量 ) CT?) (SPE) CT?) (SPE) (GLP) 2016, 67 (3): 1055-1062. (Zhao Xiaoqiang, Wang Tao. Batch process fault 
误 报 率 16.48 8.16 14.85 6.36 0 diagnosis based on TGNPE algorithm [J]. CIESC Journal, 2016, 67 (3): 
漏 报 率 18.91 15.53 14.43 12.97 8.96 1055-1062. ) 
检测 率 81.09 84.47 85.57 87.03 91.04 [4] Ben IK, Liama M, Weihs C. Variable window adaptive kernel principal 
图 4~ 7 为 采用 DPCA 和 GLNPE 对 故障 1 的 监控 ,从 图 中 component analysis for nonlinear nonstationary process monitoring [J]. 
可 以 明显 看 出 DPCA 和 GLNPE 均 有 误 报 、 漏 报 产生 。 虽 然 Computers & Industrial Engineering, 2011, 61 (3): 437-466. 
DPCA 在 一 定 程度 上 可 以 解决 动态 性 ; 但 从 表 2 可 以 看 出 [5] Chen J, Liu K C. On-line batch process monitoring using dynamic PCA 
GLNPE 的 误 报 率 、 漏 报 率 较 低 、 检 测 率 高 于 DPCA， 这 是 因为 and dynamic PLS models [J]. Chemical Engineering Science, 2002, 57 (1): 
DPCA 并 不 能 很 好 地 保持 非 线 性 数据 过 程 中 的 全 局 、 局 部 结构 ， 63-75. 
只 是 对 线性 过 程 的 动态 特性 处 理 .GLNPE 不 仅 保 持 了 全 局 结构 [6] Lee J M, Changkyoo Y, In-Beum L. Fault detection of batch processes 
言 息 , 而 且 可 以 有 效 地 提取 隐藏 在 高 维 空间 中 的 低 维 流 形 结构 ， using multi-way kernel principal component analysis [J]. Computers & 
虽然 不 是 对 非 线性 最 有 效 地 解决 途径 ， 但 通过 保持 非 线性 数据 Chemical Engineering, 2004, 28 (9): 1837-1847. 
的 全 局 、 局 部 信息 ， 故 障 检测 率 明 显 优 于 传统 直接 忽略 非 线 性 [7] Peng Kaixiang, Zhang Kai, Li Gang, ef al. Contribution rate plot for 
的 方法 。 表 3 为 故障 2 的 监控 统计 结果 ， 从 表 3 也 可 以 看 出 nonlinear quality related fault diagnosis with application to the hot strip 
GLNPE 比 DPCA 具有 更 好 的 监控 效果 ， 进 一 步 证 实 了 保持 全 mill process [J]. Control Engineering Practice, 2013, 21 (4): 360-369. 
局 -局 部 信息 对 非 线性 过 程 监控 的 有 效 性 。 [8] Fan Jicong, Wang Youqing. Fault detection and diagnosis of non-linear 
图 8、9 为 GLRGMM 对 故障 1 和 故障 2 的 监控 , 可 以 看 出 non-Gaussian dynamic processes using kernel dynamic independent 
GLRGMM 的 监控 效果 优 于 DPCA fll GLNPE. GLRGMM 可 以 component analysis [J]. Information Sciences, 2014, 259 (3): 369-379. 
检测 出 所 有 故障 1 的 变量 ， 没 有 漏 报 产生 ， 误 报 率 仅 为 0.5%， [9] Lu Ningyun, Yao Yuan, Gao Furong, et al. Two-dimensional dynamic PCA 
检测 速度 快 ; 对 故障 2 BA Dido, [EIC VG] 3 A 91.0496, for batch process monitoring [J]. AIChE Journal: Chemical Engineering 
故障 敏感 性 高 。 这 是 因为 GLRGMM 不 仅 最 大 程度 地 保持 了 非 Research and Design, 2005, 51 (12): 3300-3304. 
线性 数据 的 全 局 、 局 部 结构 ， 而 且 考 虑 了 数据 的 动态 性 ， 对 非 L0] Wang Ting, Wang Xiaogang Zhang Yingwei, er al. Fault detection of 
线性 和 动态 特性 过 程 具 有 明显 的 优越 性 。 nonlinear dynamic processes using dynamic kernel principal component 


analysis [C]// Proc of the 7th World Congress on Intelligent Control and 


5 ”结束 语 


Automation. 2008: 3009-3014. 


强 非 线性 、 动 态 性 是 间歇 过 程 的 固有 特性 ， 本 文 提 出 的 基 [11] Zhang Yingwei, Li Shuai, Teng Yongdong. Dynamic processes monitoring 
于 GLRGMM 的 监控 策略 , 不 仅 通 过 改进 NPE 对 非 线性 数据 进 using recursive kernel principal component analysis [J]. Chemical 
行 全 局 一 局 部 结构 保持 ， 而 且 通 过 引入 正则 项 对 高 斯 模型 进行 Engineering Science, 2012, 72 (16): 78-86. 
在 线 更 新 ， 最 后 整合 全 局 一 局 部 概率 指标 进行 在 线 监 控 。 通 过 [12] Jia Mingxing, Chu Fei, Wang Fuli, et al. On-line batch process monitoring 
青霉素 发 酵 过 程 的 仿真 实验 结果 表明 该 方法 比 DPCA、GLNPE using batch dynamic kernel principal component analysis [J]. 
有 更 好 的 监控 效果 ， 降 低 了 误 报 率 、 漏 报 率 ， 提 高 了 故障 监控 Chemometrics and Intelligent Laboratory Systems, 2010, 101 (2): 110-122. 
的 快速 性 和 精确 性 。 [13] Hu Kunlun, Yuan Jingqi. Multivariate statistical process control based on 


multiway locality preserving projections [J]. Journal of Process Control, 


201811.00169v1 


chinaXiv 


2008, 18 (7-8): 797-807. 

[14] Hu Kunlun, Yuan Jinggi. Statistical monitoring of fed-batch process using 
dynamic multiway neighborhood preserving embedding [J]. Chemometrics 
& Intelligent Laboratory Systems. 2008, 90 (2): 195-203. 

[15] Zhang Muguang, Ge Zhiqiang, Song Zhihuan, et al. Global-local structure 
analysis model and its application for fault detection and identification [J]. 
Industrial & Engineering Chemistry Research, 2011, 50 (11): 6837-6848. 

[16] Yu Jianbo. Local and global principal component analysis for process 
monitoring [J]. Journal of Process Control, 2012, 22 (7): 1358-1373. 

[17] Tong Chudong, Yan Xuefeng. Statistical process monitoring based on a 
multi-manifold projection algorithm [J]. Chemometrics & Intelligent 
Laboratory Systems, 2014, 130 (15): 20-28. 

[18] Yu Jie, Qin S J. Multimode process monitoring with Bayesian 


inference-based finite Gaussian mixture models [J]. AIChE Journal, 2008, 


赵 小 强 ， 等 : 一 种 基于 GLRGMM 的 间 歌 过 程 在 线 监 控 策略 7 


ChinaXiv 合 作 期 刊 


ES EE 


54 (7): 1811-1829. 

[19] Chen Tao, Zhang Jie. On-line multivariate statistical monitoring of batch 
processes using Gaussian mixture model [J]. Computers & Chemical 
Engineering, 2010, 34 (4): 500-507. 

[20] Yu Jianbo. Semiconductor manufacturing process monitoring using 
Gaussian mixture model and Bayesian method with local and nonlocal 
information [J]. IEEE Trans on Semiconduct, Manufact, 2012, 25 (3): 
480—493. 

[21] Kourti T. Multivariate dynamic data modeling for analysis and statistical 
process control of batch processes, startups and grade transitions [J]. 
Journal of Chemometrics, 2003, 17 (1): 93-109. 

[22] Birol G, Ündey C, Cinar A. A modular simulation package for fed-batch 
fermentation: and Chemical 


penicillin production [J]. Computers 


Engineering, 2002, 26 (11): 1553-1565. 


